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ABSTRACT 

The Taurus-Auriga molecular complex (TMC) is the main laboratory for the study 
of low mass star formation. The density and properties of interstellar dust are expected 
to vary across the TMC. These variations trace important processes such as dust 
nucleation or the magnetic field coupling with the cloud. In this article, we show how 
the combination of near ultraviolet (NUV) and infrared (IR) photometry can be used 
to derive the strength of the 2175 A bump and thus any enhancement in the abundance 
of small dust grains and PAHs in the dust grains size distribution. This technique is 
applied to the envelope of the TMC, mapped by the GALEX All Sky Survey (AIS). UV 
and IR photometric data have been retrieved from the GALEX-AIS and the 2MASS 
catalogues. NUV and K-band star counts have been used to identify the areas in the 
cloud envelope where the 2175 A bump is weaker than in the diffuse ISM namely, the 
low column density extensions of L1495, L1498 and L1524 in Taurus, L1545, L1548, 
L1519, L1513 in Auriga and L1482-83 in the California region. This finding agrees 
with previous results on dust evolution derived from Spitzer data and suggests that 
dust grains begin to decouple from the environmental galactic magnetic field already 
in the envelope. 
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Molecular clouds form in large structures or molecular 
complexes in specific areas of the Galaxy. The Taurus-Auriga 
Molecular Complex (TMC) covers about 20°x20° in the sky. 
The TMC is the prototype of filamentary molecular complex 
where star formation proceeds less efficiently than in mas¬ 
sive clusters, such as the L1630 cloud in the Orion molecular 
complex (Lada 19922). The large filamentary structure of 
the TMC is reminiscent of the galactic open clusters distri¬ 
bution. 


Molecular clouds are not isolated entities in the galactic in¬ 
terstellar medium. They are produced at locations where dif¬ 
fuse gas in the Galaxy is gathered due to the action of galac¬ 
tic stresses basically, spiral waves exciting Parker-Jeans in¬ 
stability modes (Elmegreen 1990, Franco et al. 2002) though 
eventually accretion from the intergalactic medium cannot 
be neglected (Mirabel 1982, Wang et al. 2004). 


1 INTRODUCTION 


tion field of the Galaxy that produces photo evaporative 
flows on the surface. The hardness of the ultraviolet (UV) 
radiation field around the TMC can be readily inferred from 
the SPEAR/FIMS imaging spectrograph observations. The 
cloud core and the halo region are clearly identified (Lee et 
al. 2006). The scattered far UV(1370A - 1670A) radiation 
is seen in, and beyond, the halo region (where A v < 1.5 mag) 
as plumes in photo-evaporative flows, on the scale of the 
SPEAR/FIMS resolution: 5-10 arcmin or 0.2-0.4 pc . Map¬ 
ping the diffuse emission and the extinction variations at 
those optical depths provides important information on the 
physics of the cloud boundary that, contrary to the cloud 
interior, is not dominated by the gravitational field but by 
more subtle stresses. Moreover, there is recent evidence of 
additional heating sources in the edge of the TMC, since the 
diffuse galactic UV radiation field cannot solely account for 
the degree of excitation of the molecular Hydrogen (Gold¬ 
smith et al. 2010). 


Molecular clouds are submitted to the ionising radia¬ 
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Cloud haloes may provide important clues on the trans¬ 
mission/reflection of long wavelength modes of the galactic 
magnetic field, as those hypothesised by Gomez de Castro 
V Pudritz (1992, hereafter GdCP92) to describe the global 
properties of the Taurus star-forming region. The geometry 
of the magnetic field in Taurus has been traced by optical 
and infrared polarisation measurements (Vrba et al. 1976, 
Moneti et al. 1984, Heyer et al. 1987, Goodman et al. 1990, 
Chapman et al. 2011). The observations indicate that the 
field is well ordered on the scale of the whole cloud, but that 
it also contains a significant disordered component on the 
smaller scales (Goodman et al. 1990). Zeeman measurements 
of the field strength in molecular clouds indicate that mag¬ 
netic and gravitational energy density are comparable. Thus, 
magnetic fields must play a fundamental role in the dynam¬ 
ics of the cloud and the formation of structures such as the 
ubiquitous cores and filaments. The orientation of the mag¬ 
netic field with respect to the dense filament substructure 
can vary from perpendicular to parallel alignments depend¬ 
ing upon the filament (Vrba et al. 1976, Heyer et al. 1987, 
Goodman et al. 1990). A recent evaluation of the field rel¬ 
evance in the gas flow suggests that magnetic fields control 
the mass flow in the TMC diffuse envelope. Heyer V Brunt 
( 2012 ) studied the alignment of the velocity anisotropy with 
the magnetic field over varying physical conditions and envi¬ 
ronments to conclude that the velocity anisotropy is aligned 
with the local, projected mean magnetic held direction in 
low surface brightness 12 CO emission areas, corresponding 
to regions of low visual extinction and presumably, low gas 
volume density. This would be consistent with the presence 
of a large scale Alfven wave channelling the dynamics of the 
TMC, as suggested by GdCP92. 

The coupling between the clouds and the ambient held 
relies on the dust grains ionisation and mass spectrum. Dust 
grains contain most of the mass of the charged particle com¬ 
ponent in the molecular clouds interior. For typical clouds, 
the cut-off wavelength for the coupling with the ambient 
held decreases from 1 pc to 0.1 pc when dust grains are con¬ 
sidered (Nakano, 1998). Also, the size of the dust grains is 
relevant for the resonance frequencies to Alfven waves prop¬ 
agation and the neutral-charge particles coupling (Pilipp et 
al. 1987). The near ultraviolet bump (the 2175 A bump) 
is the most sensitive feature to the presence of small dust 
grains (see i.e. Draine 2003) and hence, to plasma held cou¬ 
pling. Therefore, though the use of spectral tracers able 
to penetrate the dusty environment (infrared, X-ray, radio¬ 
wavelengths) is useful to study the gravitational collapse and 
the cloud fragmentation, important pieces of information 
leading to cloud formation and support are best accessed 
through UV observations. 

The GALactic Evolution eXplorer (GALEX) has sur¬ 
veyed the envelope of the TMC complex permitting a study 
of the dust spatial distribution with a sensitivity signif¬ 
icantly higher than the conventional means (molecular or 
infrared mapping). Moreover, as the TMC is a well studied 
region, comprehensive mapping in various molecular probes 
such as CO (Ungerechts V Thaddeus 1987, Heyer et al. 1987, 
Goldsmith et al. 2008, Davis et al. 2010), NH 3 (Benson V 
Myers, 1983, Gaida et al. 1984, Olano et al. 1988, Ladd et 
al. 1994) as well as in the infrared range (Froebrich et al. 
2007, hereafter F07, Lombardi et al. 2010, hereafter LLA10) 
is available allowing to study the relation between gas and 


dust density in various bands. For instance, the comparison 
between the relative extinctions in the GALEX NUV band, 
centred at the 2175 A bump, and in the K band from the 
2MASS survey (LLA10) permits to evaluate variations in 
the average dust grains size over the TMC in a very efficient 
manner. Evidence of dust grain growth in dense molecular 
gas filaments compared with the diffuse interstellar medium 
(ISM) have been recently, reported from infrared studies (see 
e.g. Ysard et al 2013 study for the LI506 filament in Taurus 
or Flagey et al. 2009 for a Spitzer based study). 

The UV bump is, by far, the strongest spectral fea¬ 
ture in the extinction curve but its source remains uncer¬ 
tain (see e.g. the review by Draine 2003). Though small 
graphite grains were proposed initially as the main source 
of the bump, the baseline today are polycyclic aromatic hy¬ 
drocarbon (PAH) molecules (Weintgarner V Draine, 2001), 
that share with a graphite sheet, a similar structure in terms 
of the distribution of the carbon atoms. PAHs are required 
to reproduce the observed infrared emission (Leger V Puget, 
1984) and they somewhat represent the extension of the dust 
grain size distribution into the molecular domain. 

In this article, we derive a method to determine the 
strength of the 2175 A bump from a combination of extinc¬ 
tion measurements in the GALEX near UV (NUV) band 
and the 2MASS K infrared band (see Sect. 2 ). In Section 3, 
we describe the GALEX All Sky Survey of the TMC, the 
quality and main characteristics of the data. In Section 4, 
we describe the method used to measure the extinction in 
the NUV band in the TMC and the results are analysed 
in Section 5. It is shown that the strength 2175 A bump 
decreases with respect to the average ISM values in the 
filaments gathering areas, namely, the low column density 
extensions of L1495, L1498 and LI524 in Taurus, LI545, 
L1548, L1519,L1513 in Auriga and L1482-83 in the Cali¬ 
fornia region. These finding agrees with previous results on 
dust evolution derived from Spitzer data. This indicates that 
within the halo of the molecular clouds, the dust grains are 
already rather large, affecting the coupling between the dust 
and the environmental magnetic field. A brief summary is 
provided at the end of the article. 


2 A METHOD TO MEASURE THE 

STRENGTH OF THE 2175 A BUMP FROM 
A nuv AND A k 

The extinction curve provides important clues about the 
properties of dust in space and the grains size distribution 
(see, for instance the review by Draine 2003 and references 
therein). The shape of the extinction curve is often modelled 
using to polynomial fits that extend from the far UV to the 
infrared using as free parameter Ry = Ay /E(B — V). In 
this work, we use the widespread extinction model derived 
by Fitzpatrick V Massa 2007 (hereafter FM07). 

The curve is well behaved from the 3 /xm to 0.3 /xm; it 
allows a simple parametrisation in the infrared that extends 
smoothly into the optical range so, 

£_ 1 + ±*(*-V, ( 1 ) 

with k(K-V) = (—0.83 + 0.63R)A -1 ' 84 —R and A = 2.2 /jm 
for the average extinction law according to FM07. Hence, 
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4^ = 1 + -h-0.19 -0.85i?) ( 2 ) 

Ay H 

This smooth trend is broken in the UV range by the 
presence of the UV bump; five coefficients are required for 
the mathematical representation of the extinction law in 
the UV (see Eq. 3), two of them to fit the UV bump to 
a Lorentzian function of (1/A). Following FM07, the extinc¬ 
tion law in the NUV band is parametrised as, 

k( A -V) = E( A - V)/E(B -V) = (3) 

Cl + c 2 x + c 3 D(x , Xo, 7 ) 

with x — 1/A (/mi) and, 

z 2 

D(x,Xo,T) = T~2 - 2\2 _l 2 2 ( 4 ) 

( X Z — Xq) Z + X Zn ( Z 

Parameters ci, c 2 , c 3 , xo and 7 are constants that depend on 
the specific line of sight. The area of the bump is given by 
dbump = 7 r 03 /( 27 ) and its maximum intensity by /bump = 
c 3 /j 2 . 

The width of the bump varies considerably from line of 
sight (LoS) to LoS with an average value of 0.993/xm -1 but 
its central wavelength remains estable (Fitzpatrick V Massa, 
1986). Absorption by small graphite particles was first hy¬ 
pothesised, as the main source of this strong absorption how¬ 
ever, variations in the graphite grain shape and size should 
produce variations both in 7 and xo (Draine V Malhotra 
1993). Since the molecular structure of PAHs molecules is 
very similar to a portion of a graphite sheet, PAHs are a 
natural extension to the graphite hypothesis and are nowa¬ 
days considered the main source of the bump (see Draine 
2003 review for more details). 

By definition, 

*^L = l + l H NVV-V) ( 5 ) 

where k(NUV — V) is weighted over the GALEX NUV band. 
Hence, 


k(NUV -V) = 



gNuv(x)(ci 


+ c 2 x + c 3 D(x , Xo , j))dx 


(6) 


— ci < gNuv > +C 2 < xgNUV > +C 3 < gNuv(x)D(x,xo,j) > 


with gNUv(x) the normalised transmittance function of the 
GALEX NUV channel Q Xl = 1680 A and x 2 = 3000 A. 
Thus, by definition, 


< gNuv >= 



gNuv(x)dx — 1 


(7) 


and, 


rx 2 

< xgNuv >— / xgNuv(x)dx — 4.49 ( 8 ) 

J X\ 

The third term in Eq. 6 depends on 7 since it influences 
D(x,x 0,7)- We have evaluated this integral for the various 
7 values in FM07; they correspond to different LoS and ex¬ 
tinction laws in the Galaxy. We have found that the integral 
can be fitted to a linear function of 7 such that, 


1 http://galexgi.gsfc.nasa.gov/docs/galex/Documents/ 


Table 1 . The standard extinction law in the GALEX NUV band 


LoS 

r(°) 

Anuv 

Ay 

Ay 

R(NUV - K)T) 

Anuv 

A k 

HD698 

3.94 

2.56 

0.09829 

9.7 

26.05 

HD3191 

2.81 

2.923 

0.07843 

7.993 

37.26 

BD+57 245 

2.97 

2.828 

0.08216 

8.156 

34.42 

BD+57 252 

2.97 

2.889 

0.08216 

8.335 

35.16 

NGC457 P34 

2.94 

2.806 

0.08149 

8.011 

34.44 

NGC457 P13 

3.11 

2.588 

0.08511 

7.784 

30.41 

NGC457 P9 

2.76 

3.089 

0.07718 

8.314 

40.03 

ISM average 

3.001 

2.735 

0.08284 

7.958 

33.01 


G) R = A V /E(B - V ) 

W R(NUV — K) = (Anuv - A k )/E(B - V ) 


< gNuv(x)D(x,x 0 , 7 ) >= (9) 

-(0.82 ± 0.02)7 + (1.38 ± 0.02) 

with rms=0.003. This tight correlation was expected given 
the functional form of D(x,x 0 , 7 ) in the GALEX band and 
the strength of the 2175 A feature. Since the area of the 
2175 A bump also depends on 7 , dbump = 71 x 3 / 27 , then, 

< gNuv(x)D(x,x 0 ,'y) >= -1.29Ab„ mp c! + 1.38 (10) 

and therefore, 

A nuv _ R + Cl + 4.49c 2 + c 3 (-1.29A bu ' mp c§ + 1.38) 

A k ~ 0.15.R — 0.19 [ ’ 

FM07 made a detailed evaluation of the extinction law 
parameters for several LoSs (see Table 4 in their article) 
and using them, we have computed the Anuv/A y, Ak/Av 
and Anuv/Ak values for these specific LoSs; also, we have 
computed them for the average ISM coefficient^] as derived 
by FM07 (see Table 1). 

In general, our R(NUV-K) is slightly larger than the 
values computed by Yuan et al. (2013), R(NUV-K) = 6.75, 
as well as their evaluation from the FM07 extinction law 
R(NUV-K) = 6.31. 

Anuv/Ak and dbump have been evaluated for all the 
LoS in FM07 (see Table 1) and plotted in Fig. 1 . Notice that 
Anuv/Ak and dbump values are well correlated unless for 
HD 3191, which seems to be a peculiar 0. For this reason, it 
has been excluded, and the resulting fit in Fig. 1 is: 

Abump = (0.106 ± 0.008) + ( 2.0 ± 0.3) ( 12 ) 

Ak 

This result is not a surprising; it just points out that 
the ISM extinction law admits a rather simple parametrisa- 
tion from the near infrared to the near ultraviolet unless for 
the strength of the bump. Therefore, for the dust grain size 
distributions observed in the ISM, the main deviations from 
this trend are caused by the UV bump. 

Likewise, the strength of the bump can be evaluated 
as a linear function of the NUV to K band extinction for 
the ISM. Note that this semi-empirical relation has been 
derived using the LoSs in FM07 that sample a broad range 


2 Parameters values for the average galactic extinction curve are: 
xo = 4.592 p m -1 , 7 = 0.922 /im -1 , c\ = —0.175, C 2 = 0.807 and 
c 3 = 2.991. 

3 HD 3191 also affects significantly the average ISM law, as 
parametrized by FM07 
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Figure 1. Relation between the area of the 2175 A bump and the 
relative NUV versus K extinction for the FM07 extinction law. 

of extinction laws (z.e., R-values and bump strengths) and 
it is independent of the physical origin of the bump be it 
PAHs, small graphite grains or any other possible source 
( e.g . fullerens, see Wada et al. 1999). 

In principle, similar expressions could be derived for 
other bands however, the availability of the near infrared 
extinction map of the TMC (LLA10) makes the Anuv/Ak 
rate, well suited for our purpose. Note that emission from 
average warm dust grains dominate the K-band while ab¬ 
sorption by PAHs or alike dominates the NUV band. 

Mid-infrared images have also been used to test dust 
evolution by comparing the emission from PAHs and very 
small grains with the radiation from large dust grains (see 
i.e. Flagey et al. 2010). However, the NUV range is ideally 
suited for this purpose given its high sensitivity to small 
particles; the estimated oscillator strength per nucleon is 
nxfx/riH — 9.3 x 10 -6 for the 2175 A bump (Draine, 1989). 
For completeness, we also note that, 

Anuv — Ak = E(B — V)(0.19 + c\ + 4.49C2 — 1.29A bu 1 mp C3 

+ 1.38c 3 + 0.85R) (13) 

and thus, ( Anuv — Ak)/E(B — V ) can be parametrised in 
terms of the strength of the bump and, as above, a simple 
scaling can be made in terms of the extinction law for various 
LoS, 

^l bump = (2.4 ± 0.5) - (0.9 ± 1.5) (14) 

though the fit is significantly less accurate (see Fig. 2). 


3 THE GALEX SURVEY OF THE TMC 

The baseline of the GALEX All Sky Survey was completed 
in 2007. GALEX AIS covers 26,000 deg 2 (~ 63% of the 
Sky) and provides broadband imaging in two UV bands. 
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Figure 2. Relation between the area of the 2175 A bump and 
the ( Anuv ~ Ak)/E(B — V) for the FM07 extinction law. 


GALEX obtained 197 images on the Taurus molecular cloud 
with a total coverage of ~ 200 deg 2 (see Fig. 3); note that 
the GALEX field of view is circular with radius 0.6°. The 
survey avoids the central region of the TMC where TMC-1 
is located, though some pointings are made around TMC- 
2 . Some molecular cores like L1495, L1498, L1544, L1515, 
LI548 and LI552 are partially covered in the survey. 

The GALEX mission provides as output products for 
each tile (or image in the AIS survey) several files contain¬ 
ing the pipeline processed images in the FUV and NUV 
bands, the intermediate calibration files and the catalogue 
of sources identified in each pointing; the sources in the cat¬ 
alogue are identified with the SExtractor procedure (Mor¬ 
rissey et al. 2007) and their FUV and NUV magnitudes are 
provided. Typically, the number of sources in the catalogue 
outnumbers by a factor of 2-3, the number of sources that 
can be identified as such from a simple inspection of the im¬ 
ages. The detection procedure clearly suffers from overesti¬ 
mation of actual sources. The number of spurious detections 
is very large. 

As the TMC is close to the galactic plane, a 
cross-identification with galactic sources from the Fourth 
USNO CCD Astrograph Catalog - U.S. Naval Observatory 
(UCAC4; Zacharias et al. 2013) and the 2MASS surveys 
(Strutskie et al. 2006) has been carried out. Typically, for 
each of our GALEX fields, there are a factor of 3 more 
sources in the 2 MASS survey and a factor of 10 in the 
UCAC4 survey; thus, they outnumber by far the UV sources 
identified by the GALEX pipeline and can be used to check 
the reliability of the GALEX identification. To complete the 
cross-identification of our sample with sources at different 
wavelengths, we have also cross-correlated our list with the 
WISE all-sky catalogue (Cutri et al. 2012) using the same 
criteria than for the other two databases. 

Magnitude completeness limits for 2MASS are 15.8, 
15.1 and 14.3 mag in J, H and K s bands, respectively. These 
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Figure 3. Fields observed by GALEX (AIS GR5) towards the 
Taurus Molecular Complex (open circles) are overlaid on the 
map of hydrogen column density ( N #) produced by LLA10 from 
2 MASS star counts. 


limits are well above the magnitude of M-type stars at the 
distance of the Taurus-Auriga molecular complex. J mag¬ 
nitude for a 4 Myr old M5 star is ~ 5.5 (e.g. Siess et al. 
2000). At a distance of 140 pc, the magnitude of such star 
J ~ 11.5 mag. Even in the case of notably ISM extinction 
Aj = 4 mag, the star would be detected by 2MASS at the 
distance of the molecular cloud. Similarly, UCAC4 is com¬ 
plete from the brightest stars to about R = 16 mag. This 
means that UCAC4 is complete for mid-M stars at the dis¬ 
tance of the Taurus cloud for extinctions Ay < 2 mag. As 
we do not expect large extinctions for our stars, because 
they were selected from the UV, we can conclude that we 
are complete for members of the star-forming region for the 
entire luminosity function, from the more massive to the less 
massive stars. Completeness limit for a typical 100s expo¬ 
sure time GALEX/AIS image in NUV and FUV is ~ 20 mag 
in both cases (Bianchi et al. 2013). 

Bona fide GALEX sources have been identified by hav¬ 
ing a 2MASS counterpart within a search radius of 3 arc- 
sec. This search radius was carefully selected after a precise 
study of the shift between 2MASS and GALEX sources in 
the TMC (see Gomez de Castro et al 2011 for details). After 
the cross-correlation, we kept a total of 163,313 UV sources 
as reliable detections. All them have been detected in the 
NUV band but only 10% of the sources have a FUV coun¬ 
terpart. Since this research is based on stellar statistics, we 
worked with the sample of confirmed NUV sources (163,313 
sources). Only 31 TTSs and and 63 candidates to TTSs are 
found in the area covered by this survey (Gomez de Cas¬ 
tro et al. 2015) thus, their impact in the overall statistics is 
minimum. 

In Fig. 4, the surface density of stars towards the TMC 
is plotted; the location of the main molecular cores as well as 
the brightest UV sources in the field (from Lee et al. 2006) 
are marked. The NUV surface density of stars is overlaid 
on the large scale extinction map of the area obtained by 



Figure 4. Density of NUV GALEX sources in the TMC. The 
densities are colour coded in stars per 3 arcmin 2 (see lateral bar). 
The 2MASS extinction map (LLA10) is overlaid for reference. The 
location of the main molecular cores and the main UV sources in 
the field (asterisks) are marked. 

Lombardi et al (2010) from the 2MASS survey of the TMC. 
A comparison between both maps (see also Fig. 4 in LLA10), 
points out that most of the area surveyed by GALEX has 
Ak ^0.2 mag that roughly correspond^] to Ay ^2.2 mag 
(or Nh ^ 4 x 10 21 atoms cm -2 ). The density of stars drops 
to zero at higher gas columns. For instance, the extension of 
some well known filaments, such as L1495 or L1498, shows in 
the stellar density map because the density of stars drops to 
0. Also, some diffuse structures are noticeable. For instance, 
at 5 ga i = —13°, the stellar density varies from 5-6 stars per 

3 arcmin 2 at Z ga i = 165° (westwards of the TMC) to 2-3 stars 
per 3 arcmin 2 at Z ga i = 180°; this feature passed unattended 
in previous CO maps of the area. 

4 METHODS 

4.1 The Wolf method 

The Wolf method is designed to measure extinction using as 
reference for the stellar luminosity function a nearby stellar 
field, assumed to be unextincted. The luminosity function 
is defined in terms of the cumulative star counts distribu¬ 
tion per apparent magnitude bin, in the fiducial field (Wolf 
1923: see also Gorbikov V Brosch, 2011 for a recent appli¬ 
cation of the method to study the extinction law towards 
the north celestial cap). The method is based on the fact 
that all the stars behind the same absorbing cloud are ex- 
tincted by the same amount, A\. As a result, the cumulative 
star counts distribution as a function of the apparent magni¬ 
tude, N(m( A)), is homogeneously shifted to lower apparent 

4 The average ISM extinction law, as described in Sect. 2, has 
been used for this scaling. 
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Figure 5. Determination of the extinction using the Wolf method 
for tile AIS190_sg27. The plot illustrates the cumulative counts 
distribution in AIS190_sg27 (solid line) compared with AIS189_28, 
the reference tile (dashed). 


magnitudes in any dusty area. The shift of the cumulative 
counts distribution is a direct measure of the relative ex¬ 
tinction between the fiducial field and any given field. The 
process is illustrated in Fig. 5 for tile AIS190_sg57. A control 
tile, tile AIS189_28, has been used for calibration purposes 
since extinction is not apparent neither in the UV data (see 
Fig. 3) nor in the infrared data (see LLA10). The NUV cu¬ 
mulative distribution in AIS190_sg297 is clearly shifted with 
respect to the reference field pointing out a NUV extinction 
of Anuv = 2.8 mag. 

Following this procedure both Anuv and Ak values 
have been computed however, the Anuv/Ak ratios obtained 
are unusually small and inconsistent with the interstellar 
extinction law ( e.g . FM07). These small ratios are obtained 
because the GALEX AIS fields contain areas with and with¬ 
out dust clouds; this geometric dilution together with the 
clumpiness of the ISM produces a weighted average of the 
A^uv/Ak value over the GALEX AIS field of view (see Ap¬ 
pendix A for more details). Therefore, it is critical to work 
with spatial resolutions that adapt well to the characteristic 
scales of the ISM structures to be studied such as filaments, 
globules etc, typically about 0.1-0.4 pc. 

Unfortunately, the Wolf method requires a good statis¬ 
tics to sample properly the luminosity function forcing the 
resolution element to be much larger than 0.1-0.4 pc (2.5-10 
arcmin scales at the TMC distance of 140 pc). A fine tun¬ 
ing between resolution and statistics is required. The star 
counts method is better suited for this purpose. 


4.2 Star counts 

Since the seminal work by van Rhijn (1929), star counts 
have been used to measure extinction in galactic fields. A 
detailed description of the method is given by Bok V Cord- 
well (1973). This technique has been extensively used in as¬ 


tronomy and has been successfully applied to measure the 
gas column to Taurus using the R plate of the Palomar Ob¬ 
servatory Sky Survey (Cernicharo V Bachiller, 1984). 

Basically, the apparent magnitude of a given star in the 
NUV band, m(NUV ), is given by: 

m(NUV) - M(NUV) = -5 + 5 log — + ,4nuv (15) 

pc 

where M(NUV) is the absolute magnitude, d is the distance 
and Anuv is the extinction in the NUV band caused by the 
intervening gas. In its most usual application, star counts 
in a given field are compared with the predictions for an 
unextincted, nearby field and the extinction is given by, 

A\ = \og(Nx/N\)/b\ (16) 

where N\ are the observed counts, are the ex¬ 

pected counts from a non-extinguished field, and b\ — 
d\og N(m\)/dm\ (with m\ the apparent magnitude at 
wavelength A) is a measure of the slope of the luminosity 
function in the area under study. 

The NUV luminosity function has not been determined 
for the Galaxy thus we have derived it in this work for the 
TMC region. Firstly, we have determined it for the reference 
field (AIS 189_28), which we consider unextincted and has 
2,340 sources. We have calculated it to be well fitted by, 


log N(rriNuv) — (0.294 d= 0.002 )rriNuv — 2.55 A 0.04 


with RMS = 0.046 and thus, bNuv = 0.294 (see Fig. 6). 
We also have tested whether there are significant variations 
over the TMC field. As shown in the bottom panel of Fig. 6, 
the slope remains constant (within the error bars) pointing 
out that only minor variations are detected over the field, if 
any. 

The accuracy of the Anuv values derived from the star 
counts method is (see Dickman 1978), 


5NUV _ Nnuv_ x - i / Nnuv_ A N^ uv 1/2 
NUV ~ [ ° g N NUV ) 1 N nuv N^ uv J 

with, 


SAnuv— 


1 

bNUV 


( 


Nnuv + N^uy 
NnuvN nuv 


xl/2 


(18) 


Thus, for Anuv — 0 (or N^ uv — Nnuv ), SAnuv — 
0.1 mag. 

The maximum Anuv is observed in tile AIS190-43 with 
Nnuv — 136, and thus Anuv — 3.2 mag with 5Anuv 
0.8 mag. 

If the statistics are good then both the Wolf and 
star-counts should produce very similar results. Av¬ 
erage Anuvs have been calculated by both methods for all 
the GALEX fields and, as shown in Fig. 7, the only dis¬ 
crepancy is at high extinctions (low star counts) since the 
Wolf method produces slightly higher extinctions than the 
star-counts method. The regression line is: 


Anuv (Wolf) = (1.10 A 0.03) Anuv (Counts) + (0.08 A 0.04) 


with rms — 0.27. 
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Figure 7. Comparison between the A nuv values obtained by 
direct star counts and by using the Wolf method. The straighline 
marking A nuv ~ Wolf method equal to A nuv -Star counts is 
plotted for comparison. 


being Anuv, the extinction determined with the star counts 
method, Eq. 16, and A^ uv , the final value of extinction. To 
evaluate this correction, we have made use of the Besancon 
model of the Galaxy (Robin et al. 2003; Czekaj et al. 2014, 
see also URL: model.obs-besancon.fr), and the Ak map of 
the TMC evaluated by LLA10 (see details in Appendix B). 
The correction factor is found to be, 

dA m (0.0640 ± 0.0005)6^ (°) - (0.38 ± 0.01) (20) 

(with rms — 0.025). 

The variation of the stellar luminosity function between 
Zgai = 165° and l ga \ = 180° is negligible according to the 
Becanson standard model. 


Figure 6 . Stellar luminosity function in the TMC in the NUV 
band. Top: NUV luminosity function in the reference field, bot¬ 
tom: luminosity function derived from all NUV sources in the 
AIS TMC survey. The best fit to the slope is shown by a dashed 
line on the fitted region. 


4.3 Galactic latitude correction 

Stellar density declines significantly with galactic latitudjfl 
To correct for this effect, we have computed the variation of 
N^uv with galactic latitude and derive a correction factor, 
dA, such that, 

4nuv = Anuv — gL4nuv (19) 


5 As pointed out in the text, we have not found significant vari¬ 
ations in the luminosity function, as a simple inspection of Fig. 6 
shows. 


5 RESULTS: VARIATION OF A NU v/A K OVER 

THE TMC 

Cores and filaments in molecular clouds have typical physi¬ 
cal scales of 0.1-0.4 pc that correspond to 2.5-10 arcmin at 
the TMC distance (140 pc). Unfortunately, the UV statistics 
is too poor at resolutions of 3 arcmin 2 (see Fig. 4). Thus, 
we have followed an incremental approach using two inter¬ 
mediate bin sizes of 24 arcmin 2 and 12 arcmin 2 . 

NUV extinctions have been calculated by applying the 
star counts method (see Sect. 4) including the galactic lat¬ 
itude correction derived in Appendix B. To compute K- 
extinctions we have resourced directly to Froebrich et al. 
(2007) map. The map has been rebinned and re-scaled to 
K-magnitude; note that F07 provide the electronic version 
of the map (FITS format) in visual extinctions (Ay), but 
these Ay values are obtained by the authors after applying a 
simple scaling from infrared to visual magnitudes assuming 
that the standard ISM extinction law holds (Draine 1989). 

As shown in Fig. 8, there is a significant variation of 
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Anuv / Ak across the region. Anuv / Ak — 33 in the dif¬ 
fuse ISM, around the TMC in agreement with FM07. How¬ 
ever, westwards, at the tail of the cloud, where the main 
filaments reside, and the column density of molecular gas 
increases, Anuv/Ak decreases pointing out the weakening 
of the UV bump and, in general, a larger average size of the 
dust grains when compared with the diffuse ISM. This trend 
is confirmed at both intermediate resolutions. 

Some of the areas mapped by GALEX have very low 
and thus, uncertain, Ak values (LLA10 estimates cL4k = 
±0.04 mag). This results in abnormally high Anuv/Ak 
rates, especially at high galactic latitudes (see Fig. 8 ). 

Experimentally, the UV bump is known to be produced 
by mixtures of sufficiently large, neutral PAHs (Steglich et 
al. 2011). Moreover, the UV-visible measurements provide 
a reliable fingerprint of the presence of specific PAHs since 
the electronic transitions of these species are detected in 
this range while the near infrared vibrational bands are not 
molecule specific. 

A recent UV-visible (3050-3850 A ) study based on 
the observations towards five lines of sight with moder¬ 
ate reddenings measured PAHs abundances two orders of 
magnitude smaller than those in the diffuse ISM (Gredel 
et al. 2011 ); the measurements were made over column 
densities of Nh — 5.7 — 9.2 X 10 21 cm -2 , slightly above 
the dominant column density of our study in the TMC, 
Nh ^ 4 x 10 21 atoms cm -2 . This study was sensitive to 
molecules such as anthracene, pyrene or benzofluorene. Lab¬ 
oratory experiments suggest that the UV bump is produced 
by larger PAHs with sizes above 50-60 carbon molecules 
(Steglich et al., 2010) thus, our results suggests that large 
PAHs are neither abundant in traslucent clouds. 

Anuv/Ak is small close to the filamentary networks 
(compare Fig 3 with Fig 8 ) where also mid infrared data 
have detected the formation of large dust grains (Flagey et 
al. 2010 ), as high density of molecular gas makes collisional 
processes more efficient assisting grain growth. In addition, 
the large density of dust in the filaments effectively shields 
the growing mantles against the environmental UV radia¬ 
tion. A detailed comparison between Figure 4 and 8 , shows 
that the bump is less prominent towards the tail of the fil¬ 
aments leading to L1495 and L1498. The first results from 
the Herschel survey of the Gould Belt in Taurus, show a 
fraction of this area: Nl-3 for the filament leading to L1495 
and Sl-4 for that leading to L1498. The analysis of Herschel 
data shows that these filaments consist really in a hierarchy 
of cores of molecular gas with sizes in the range 0.0242.7pc 
that covers the full core mass function from unbound clumps 
to (gravitationally bound) pre-stellar cores, being all part of 
the same population (Kirk et al. 2013). Our detection of 
increased average size of dust grains around this area is con¬ 
sistent with the enhancement of the molecular gas and dust 
density that drives the formation of molecular cores. 

In the global scenario developed by GdCP92, matter 
in the TMC is being gathered by the local action of the 
Parker-Jeans instability. Matter lifted to the halo by the ac¬ 
tion of the instability, rains back along the magnetic field 
lines reaching velocities comparable to the Alfven speed in 
the undisturbed gas layer in the galactic plane. These ve¬ 
locities agree well with the clouds and pre-main sequence 
stars velocities that roughly run eastwards (see GdCP92). 
Fresh material is thus expected to be gathered at the rear 



Figure 8. Anuv/Ak on the TMC computed for bins of 
24 arcmin 2 and 12 arcmin 2 , top and bottom panels respec¬ 
tively. The gray scale runs from Anuv / Ak = 0 (black) to 
Anuv/A y ^ 33 (white). Anuv / Ak ^ 33 corresponds to the 
predictions of the standard ISM extinction law (see Table 1). 


of the cloud, where it is, indeed. As the density increases 
with respect to the average ISM values, so does the opac¬ 
ity protecting dust grains from the ambient UV radiation 
field and favouring grain growth. In this sense, our results 
are consistent with the expectations. Unfortunately, the long 
wavelength Alfven wave mapped by Moneti et al. (1984) was 
not properly surveyed by GALEX. 


6 CONCLUSIONS 

In this work, we have shown the potentials of wide field UV 
imaging to study the distribution and properties of diffuse 
dust clouds in space. Our main contributions to the analysis 
of the UV data have been: 

• to show that the wide field UV images obtained by 
means of the GALEX AIS are able to detect very low 
dust columns that are not detected through infrared surveys 
(2MASS). These images are obtained with a rather modest 
telescope (40 cm) in 150 s (barely 2.5 minutes), in average. 

• the stars luminosity function has been derived for the 
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GALEX NUV band in the galactic anticenter region. bNuv 
has been computed. This value is suitable to be compared 
with measurements across the Galaxy. 

• a generic method to determine the relative abundance 
of small dust grains from star counts and surveys has been 
described. 

These developments have been used to detect variations 
in the extinction law across the TMC. The extinction law 
tends to be greyer in the rear of the cloud and in the net¬ 
works of filaments, than on the clean ISM. This is consistent 
with the generic picture of large dust grains being formed in 
molecular clouds and small dust grains being more abundant 
in the ISM. It is worth remarking that the areas where dust 
is being gathered are at the rear of the cloud, as expected 
if star formation at the head of the TMC is associated with 
a local storage of mass being gathered by the action of the 
Parker-Jeans instability. 

Unfortunately, the depth of the GALEX AIS is not 
enough to have good statistics in the areas around the fila¬ 
ments so, we cannot test effects like grain growth. Note that 
according to, 

A 1 nuv — l°g(A’ivc/v/^^L^v)/0-29 

and setting N l ^ v — 1, A l ^ v depends on the number of 
sources in the test field (which is the maximum) and thus, 
on the resolution. To reach spatial resolutions of 0.05 pc, rea¬ 
sonable to study the filaments structure, and A l ^ v — 4, it is 
required that N^ uv = 16 per 1.2 arcmin 2 , which represents 
a factor of 14 increase in stellar density (stars/arcmin 2 ) with 
respect to GALEX AIS. The limiting magnitude should rise 
from NUV = 22.3 mag (see Gomez de Castro et al. 2015) 
to NUV = 23.5 mag. 
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APPENDIX A: ESTIMATE OF THE IMPACT 
ON THE Anuv/Ak RATE OF THE FILLING 
FACTOR OF THE DUST CLOUDS WITHIN 
THE GALEX AIS FIELD OF VIEW 


Let Anuv and Ak be the extinctions evaluated from the 
star counts and affected by the filling factor of the extincted 
over un-extincted areas in the held of view. Let A* NUV and 
A f K the true, average values, taking into account only the 
extincted areas. Then, 

, Nnuv,' 


log 


log 


Nnuv 
N k ,max 7 A f 

^r = bKA « 


— bNuvAL , 


(Al) 


(A2) 


and, 

1 f 

6A N uv — j -logx — A J NU v ~ Anuv 

onuv 


SA k = T^- log x = A f K - A k 

ok 


(A3) 

(A4) 


with 5 Anuv and SAk the excess extinctions caused by the 
geometric filling factor of dusty areas in the GALEX held of 
view, x- With this definition, 


Anuv _ bx bNuvA f NUV — logx 
Ak bxuv b K A f K - log x 

and after some algebra, 
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making use of the bx and bNuv values determined from 
the luminosity functions (see Sect. 4.1) and assuming 
Anuv/Ak — 33.01, the average ISM value, we derive, 
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(All) 


Finally, note that ( Anuv — Ak) is basically insensitive 
to the filling factor since, 


Anuv — Ak = 
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(A13) 

(A14) 


assuming bNuv — bx- Unfortunately, this colour has a softer 
dependence on the strength of the 2175 A bump than the 
rate Anuv/Ak as shown in Sect. 2. 

This short exercise provides a first order estimate and 
neglects the effects caused by the cloud clumpiness (see, for 
instance, Lombardi 2009). 


APPENDIX B: DERIVATION OF THE 
CORRECTION FOR GALACTIC LATITUDE 

The GALEX AIS survey spans over a wide range of galactic 
latitudes (from b ga i — —5° to b ga i — —25°) and it is ex¬ 
pected that the stellar density decreases significantly over 
the surveyed area. As a result, in the evaluation of the ex¬ 
tinction: 


A\ = log(Nx/N\)/b\ 


Nx is a function of b ga i . In our work, we have used only a 
reference held, AIS189.28, located at b ga i = —6° hence it is 
necessary to determine the correction term, SAi arn bdai such 
that, 


dAx{b g ai ) 


blambda N(b ga [) 


To derive it, we have used the infrared data and a simplified 
model of the galactic stellar density distribution. 

LLA10 reported that the K-band star counts followed 
the galactic model. Both FM07 and Lombardi et al (2010) 
designed procedures to correct for this gradient so they pro¬ 
vide a reliable baseline to compare our correction with. Note 
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that no significant variations are expected in b\. Neither the 
stellar density varies significantly with galactic longitude in 
the anticenter region. 

The impact of this correction in extinction determina¬ 
tions is illustrated in Fig. Bl. Rough Ak (without SAk cor¬ 
rection ) have been computed by applying the star counts 
method on the sample constituted by all bona fidae 2MASS 
sources in any given of the GALEX fields to the depth of 
the 2MASS survey (see Fig. Bl, top panel), as well as to 
a reduced sample including only the 2MASS sources with 
NUV counterpart. Both panels display a clear gradient with 
galactic latitude, more pronounced in the top panel. At the 
bottom, F07’s infrared extinction map is shown for compar¬ 
ison (the map has been rebinned to the GALEX tiles size 
and distribution). 

According to the standard galactic model, the stellar 
density, p *, in the Galaxy can be parametrized in terms of 
a thin, p*, and a thick disk, p£, as follows: 

P* = P*a + Pb (Bl) 


with, 


R z 

Pa,b = P0,(a,6)exp(- 7 -)exp(-——) 

Ta,b Jla,b 


(B2) 


with, R and z the galactocentric cylindric coordinates in the 
disk (cylindric symmetry is assumed) and pj,a and Po,6 the 
volumetric density of stars in the center of the thin and thick 
disk, respectively. According to Ojha (2001), the parameters 
for the model are: L a = 2.8 d= 0.3 kpc, L& — 3.71 q'J kpc, 
H a — 250 pc (Kuijken &, Gilmore 1989, Haywood 1997), 
Hb = 860 =b 200 pc and Po,fe/Po,a — 0.035 =b 0.02. With these 
provisions, the total number of stars with absolute magni¬ 
tude, M, observed within a solid angle, up to a given 
depth rd, is given by, 


l 


N(m) — / &(M)p*Qr 2 dr 
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with <f>(M) the galactic distribution of stars with absolute 
magnitude, M. Let us assume that <t>(M) can be consider 
constant and equal to <l>o for all the sources in our sur¬ 
vey. This would imply that all sources follow the simple 
parametrization in Eq. (B2), irrespective of their nature, 
on scales of 1.4 deg (the resolution of the survey). With this 
approach, 
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Since the galactic longitude of our region of interest is 
lgal — 180°, z ~ (r + 7*0) tan^az), with r© the galacto¬ 
centric distance of the Local Standard of Rest (LSR) and 
(■ r,lgai,bgai ) the LSR-centric galactic coordinates. Then, 
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Figure Bl. K-band color-coded average extinction per 
GALEX AIS tile in the TMC. Top: <5Ak derived from all 
2MASS sources per GALEX tile, Middle: <5Ak derived from stars 
with NUV counterpart, Bottom: <5Ak derived from binning F07 
2 MASS based extinction map. 


Then, 
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Figure B2. Effect of the variation of stellar density with galactic 
latitude in the extinctions computed from star counts with a given 
reference field at b ga i = —6° 
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that depends on galactic latitude and depth of the survey. 
Let us denote, 
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and, 
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dA\(bgai,r d ) 


<d^ log (++ 


(BIO) 


the effect of the stellar density gradient in the extinction de¬ 
rived from star counts. Note that Eq.(BlO) is parametrized 
in terms of the reference field AIS189_28 (see Sect. 5), at 
bgai = - 6 °. 

In Fig. B2, dA\(bgai,rd) is plotted for a range of galac¬ 
tic latitudes and possible depths of the surveys according to 
Eq. (BIO). To determine the depth of the survey (and the 
proper correction to be applied), we have used three tiles, 
at different galactic latitudes: the reference tile AIS189.28 
(bgai ^ -6°), AIS 190_42 (b gai ~ -12°) and AIS169.28 
(bgai — —24°). We have cross-checked that Ak = 0 for all 
of them in the LLA10 extinction map. 

We have used Ak measurements (see Fig. Bl) to evalu¬ 
ate gL4a(—12 °, r d ) = A X (AIS190A2) - A X (AIS189.28) and 
dA x (-24°,r d ) = A X (AIS169J28) - A X (AIS 189.28). They 
are represented by dashed lines in Fig. B2 as, dA x (— 24°) 
and dAx(—12 °). 

As it is noticeable, the dashed lines meet the galactic 
model at r d ~ 800 pc, for both galactic latitudes. Thus, we 
assume that this is the depth of the GALEX AIS survey in 
this area of the Galaxy. From that, we compute the correc¬ 
tion factor to be applied to the measurements in Table 2 to 
be, 


dA = (0.0640 ± 0.0005) \b ga i(°)\ - (0.39 ± 0.01) 



185 180 175 170 155 


Figure B3. Average extinction per GALEX tile in the NUV 
band determined from star counts. Top: A^uv without applying 
the SAkuv correction, bottom: A^uv after applying the correc¬ 
tion. The location of the obscuring clouds from F07 is outlined. 
Coordinates as in Fig 1. 


(with rms = 0.025). The extinction map in the NUV band, 
after applying this correction, is displayed in Fig. B3, for 
the record. The location of the obscuring clouds, the high 
latitude cirrus and the filaments are readily shown. 
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